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Abstract 

In recent years several planets have been discovered at wide orbits (> 100 AU) around their host 
stars. Theoretical studies encounter difficulties in explaining their formation and origin. Here we 
propose a novel scenario for the production of planetary systems at such orbits, through the dynamical 
recapture of free floating planets (FFPs) in dispersing stellar clusters and stellar associations. This 
process is a natural extension of the recently suggested scenario for the formation of wide stellar 
binaries. We use N-body simulations of dispersing clusters with 10 — 1000 stars and comparable 
numbers of FFPs to study this process. We find that planets are captured into wide orbits in the typical 
range ^ few x 100 — 10^ AU, and have a wide range of eccentricities (thermal distribution). Typically, 
3 — 6 X (/ffp/1) % of all stars capture a planetary companion with such properties (where /ffp is the 
number of FFP per star in the birth clusters). The planetary capture efficiency is comparable to that 
of capture-formed stellar-binaries, and shows a similar dependence on the cluster size and structure. 
It is almost independent of the specific planetary mass; planets as well as sub-stellar companions of 
any mass can be captured. The capture efficiency decreases with increasing cluster size, and for a 
given cluster size the it increases with the host /primary mass. We also find that more than one planet 
can be captured around the same host through independent consecutive captures; similarly planets 
can be captured into binary systems, both in circumstellar and circumbinary orbits. We also expect 
planets to be captured into pre-existing planetary (and protoplanetary systems) as well as into orbits 
around black holes and massive white dwarfs, if these formed early enough before the cluster dispersal. 
In particular, stellar black holes have a high capture efficiency (> 50% and 5 — 10 x (/ffp/1) % for 
capture of stars and planetary companions, respectively) due to their large mass. Finally, although 
rare, two FFPs or brown dwarfs can become bound and form a FFP-binary system with no stellar 
host. 



1. INTRODUCTION 

Exo-planetary systems are known to exist in a wide 
variety of configurations, raising many challenges to our 
understanding of the complex processes of planet forma- 
tion. Some of these planetary systems were found to 
have massive planets orbiting their host stars at w ide or- 
bits (hundreds of AU: iLa frcnicrc ct al. 2008; Billc r et al.l 
I20TI [Ireland et al.l[20TTl: fLafrenike et al.H201lD . Such 
systems are not easily produced in current models of 
planet formation, and their origin is stil l debated (e.g. 
Dodson-Robinson et al.l 120091 : iBossI 120111 : iKratter et all 
2010; Vcras ct al. 2009). Microlensing observations find- 
ings show the possible existence of a large population 
of apparently free fioating planets (FFPs, l-Stp g per 
star; iSumi et all 120 lit ), i.e. they appear not to have 
any host staiQ, consistent with theoretical expectations 
of the production o f such ru naway planets throug h vari- 
ous processes (e.g.-lKroupa fc Bouvicr 2003; Vcra s et al.l 
[2001 iBeauge fc Nesvornvl 12011: Parker fc Quanzll2012L 
and references therein). Here we suggest a novel sce- 
nario connecting FFPs with wide-orbits planets, namely, 

|hperets@physics.technion.ac.il;thijskouwenhoven<0)gmail.com 
^ We do caution that the actual fraction of FFPs is unknown; 
it is possible that the apparently host-less planets observed have 
companions too faint to be resolved or too far away. We therefore 
scale our results throughout this paper by the FFP fraction per 
star. 



the production of wide-orbit planetary systems from the 
re-capture of FFPs in dispersing clusters. 

Recently, several studies have quantitatively shown 
that wide stellar binaries can form following the 
disper sal of a stellar cluster or stellar associa- 
tion dKouwenhoyen et al.l I2010L : hereafter paper I; 
IMoeckel fc Bate! 120101 : IMoeckel fc ClaTk^ [MTh . and a 
similar suggestion was raised for the formation of Oort 
clouds (Leyison et al. 2010). Since the ma jority of stars 
are t hought to be born in star clusters ()Lada fc Ladal 
I2003D . the majority of planetary systems should also 
form and evolve in such environments. Indeed, many 
studies explored the evolution of planets in stellar 
clusters (prior to the clusters dispersal), showing they 
could be dynamically excited by encounter with other 
cluster stars or even be expelled from their host star to 
become free floating planets (e.g. i Laughlin fc Adami 
19981: IBonnell et all [200TI: ISiSith fc Bonnell ,20011 



Adams et al.l 120061: iFreeeau et all 120061: iSpurzem et al 



20091: iMalmberget al.ri2011l : iParker fc Ouand 12012 
and references therein). Studies of the the kinematic 
properties of such FFPs suggest th ey have similar 
veloc ity dispersion as the cluster stars ([Parker fc Ouan3 
[20T1 . The suggested scenario for the formation of wide 
binaries following the cluster dispersal could therefore 
similarly produce wide planetary systems through the 
same capture process, where we assume that the FFP 
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populations is formed prior to the cluster dispersal, 
which is generally consistent with the suggested origins 
of such FFPs. 

The capture processes can be described as follows (see 
also paper I); in an evolving star cluster, stars and FFPs 
may become bound to each other as the cluster expands, 
i.e. if the gravitational influence of the other cluster 
members decreases. In order to form a binary pair in this 
way, (i) they need to be sufficiently close together, (ii) 
a star and a planet (or even two planets/brown dwarfs) 
need to have a sufficiently small velocity difference and 
(iii) the newly formed binary should not be destroyed 
by gravitational interaction with the remaining cluster 
stars, field stars or the Galactic tide. In the following we 
present the results of detailed N-body simulations which 
show that wide-orbit planetary systems around both sin- 
gle and binary stars can form at large numbers under 
realistic conditions. In addition we note that this mech- 
anism can produce FFP-binary systems, brown-dwarf 
planetary systems such as be en recently observed (e.g. 
iJavawardhana fc Ivanovl [2006D as well as planetary sys- 
tems around compact objects (white dwarfs and black 
holes). 

The paper is organized as follows. In Section §2 we pro- 
vide a brief overview of our technique and assumptions, 
which are described in details in in Paper I. In Section §3 
we describe the results of our simulations, and in Section 
§4 we present and discuss our results. We summarize our 
conclusions in Section §5. 

2. SIMULATIONS 

Generally, our models follow the same scheme as 
described in detail in Paper I; in the following we 
briefly review the simulations. Following Paper I, we 
simulate star clusters usin g the STARLAB package 
(jPortegies Zwart et al.l [2001|) . We draw N single stars 
from the IKroupal (I200lir mass function, /m(M), in the 
mass range 0.1 < M < 7M0. The lower-limit corre- 
sponds to the hydrogen-burning limit. The upper limit 
is set to the mass of the most massive star which has not 
yet evolved due to stellar evolution in 50 Myrs, the typi- 
cal timescale at the end of our simulations. We also run 
simulations with higher cut-off (20 and 50 Mq) to probe 
capture by massive stars in fast dispersing clusters. 

We perform simulations with varying N, ranging from 
small stellar systems (or sub-clumps) with = 10 to 
intermediate systems N = 20, 50, 100, 200, 500 and up 
to open cluster-sized systems {N — 1000). The cluster 
size which encloses all stars initially is taken to be related 
to the number of stars through R ~ O.IN-^^^ pc, which 
corresponds to an identical initial stellar density for all 
the simulated clusters. 

We study two sets of dynamical models to describe the 
cluster structure: Plummer models and sub-structured 
(fractal) models; these are described in detail in Paper I. 
Fractal (clumpy) models are characterized by their frac- 
tal dime nsion a, and are stru ctured following the method 
used by lAllison et al.l (|2009( ): we use two models, fractal 
clumpy clusters with a = 1.6 and homogeneous systems 
with a = 3; in this study we explore only initially ex- 
panding clusters with Q = 3/2, where Q is the virial 
ratio, defined by Q = —Ek/Ep , where Ek and Ep are 
the kinetic and potential energy of the cluster (see details 
in Paper I). 



All simulations are performed until the clusters are 
completely dissolved, which is typically of the order of 50 
Myr (we follow the largest clusters up to 200 Myrs to en- 
sure they completely dissolve), the timescale at which the 
ma jority of low-mass star clusters are destroyed (s ee, e.g. 
Tutuko v 1978; .B outloukos & Laniers 2003; B astian et al.l 
i2005 : IFalletaniloOSr and numerous others). Note that 
the lifetime of massive stars with M > 7 Mq could be 
shorter than the timescale for the cluster dispersal, and 
therefore black holes (BHs) and massive white dwarfs 
(WDs) may already form in the cluster during its disper- 
sal. Neutron stars may also form in the cluster but are 
not likely to be retained in the cluster due to their natal 
kicks. 

Binary systems and multiple systems (where the 
masses can be stellar or planetary) are identified re- 
cursively. We impose that the systems must be stable. 
The internal stability is partially guaranteed by the fact 
that they survive 50 Myr in the simula tions. Further- 
more, we impose the stability criterion bv lValtonen et all 
((2008). In our simulations we account for the tidal 
fiel d of the galaxy at the Solar neighborhood (follow- 
ing iKrouE^eT^ HqOII; we adopt a distance of 8.5 kpc 
from the Galactic c enter and take Mgai = 5 x 10^° Mq). 
iJiang fc Tremaini (|2010D show that binaries can survive 
in the Galactic tidal field if their separations are smaller 
than the Jacoby radius 
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where G is the gravitational constant, the angular ve- 
locity of the Galaxy at the solar circle, A is Oort's A 
constant, and mi and TO2 are the masses of the com- 
ponents. The value of 1.7 pc is valid for the canonical 
values for the Galactic constants and individual masses 
of 1 Mq. We do not follow the binaries longer than the 
simulation time, and therefore some of the widest sys- 
tems formed in the simulation that we consider as bound 
would be slowly disrupted in the Galactic field, later on. 
Nevertheless, 90 % of all the capture-formed systems 
we find are within their Jacoby radius, consistent with 
the expectations. Note, that wide binaries are also likely 
to be destroyed through encounters with other field stars 
or giant molecular clouds; these are not included in our 
simulations. 

The most important difference from previous simu- 
lations is the addition of a population of FFPs. We 
simulated FFP populations with one planet per star, 
/ffp = Nppp/Ni, — 1, and a more limited number of 
simulations with /ffp — 0.5 and 2. For simplicity we al- 
ways make use of single mass FFPs, taken to be Jupiter 
mass; mpFP = Mj. Since the binding potential of a 
planet and a star are completely dominated by the stel- 
lar mass, the specific choice of the planetary mass is typ- 
ically not important for their capture around stars, as 
long as TOp <^ Mi, where rrXp and are the typical 
masses of the FFPs and the cluster stars, respectively. 
Nevertheless, we also ran simulations with lower mass, 
50 Earth mass planets, mppp — 5OM0, to verify this. 

The parameters of the various cluster models are 
summarized in Table 1. Throughout the paper we 
refer to models by their parameters, e.g. model 
F_N20_P0.5_M0.16 corresponds to a fractal cluster of 20 
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TABLE 1 

Parameters of the cluster simulation runs used in our 
analysis 



Parameter 


Value 


Structure 


Fractal (F; a = 1.6), 
Homogeneous (H; a = 3) 


Number of stars (A^^) 


10, 20, 50, 
100, 200, 500, 1000 


Planets to stars ratio (/ppp) 


0.5, 1, 2 


Mass Function 


Kroupa MF, 
truncated at m = 7 Mq 


Planetary mass 


Mj,50Me 


Cluster Radius (R)) 


o.iA^y^ pc 


Simulation time 


50-200 Myrs 



Stars and 10 FFPs (/ffp = 0.5) of mass 0.16 Mj; if not 
explicitly written, /ffp ~ 1 and AfpFp = Mj. 

3. RESULTS 

In the following we describe the various populations 
of planetary systems we obtain for the different cluster 
models we explore. Table 2 summarizes the results for 
the total fractions of multiple systems formed in our sim- 
ulations. 

3.1. Fractions of capture-formed systems 

As can be seen in Fig. 1, the fractions of capture- 
formed stellar-binaries and planetary systems follow sim- 
ilar dependence on cluster sizes and cluster structure (ho- 
mogeneous vs. fractal). The binary fraction in homoge- 
neous clusters is inversely proportional to the number 
of stars in the cluster, but shows a weaker dependence 
(/cap oc iV~°-^^) on the cluster size than suggested by 
the theoretical expectation of /cap ~ from the the - 

ory of soft binaries freeze-out (jMoeckel fc Clarkil201lD . 
This model predicts a single binary to form per cluster, 
independent of the number of stars in the cluster. The 
fractal clusters also show a decrease in capture efficiency 
with increasing cluster size, with a stronger dependence, 

/capOciV-O-35. 

We can now integrate the observed fractions (from 
our simulations) together with the observed cluster mass 
function dN/dN^ oc A^r^ dLada fc Ladal [20031) in the 
range 10 < iV^, < 1000, and calculate an overall ex- 
pected fraction of capture-formed systems. For the ho- 
mogeneous clusters we find fractions of f^^ ~ 0.069 and 
/planet ^ 0.058 for the stellar and planetary systems, 
respectively; similarly, for the fractal clusters, we get 
/bT„P~ 0.055 and /p^P,,^ 0.036. _ 

The observed wide binary fraction (with semi-major 
axis SMA, a > 10'^ AU) is ~ 0.02 fRaghava n et al.l 
120101: IMoeckel fc Clarkd I201ll ): recent observations fo- 
cusing on binaries at larger sep arations suggest that the 
fraction could be even higher (IShava fc Olling|[20Tll) . if 
these are accounted for. In the simulations a large frac- 
tion (~ 80 % and 50 % for the homogeneous and 
fractal clusters, respectively) of the capture-formed mul- 
tiple stellar systems reside in this range. We therefore 
conclude that the observed fraction of wide binaries in 
our simulations is generally consistent with observations. 
Taken at face value, we may therefore expect the effi- 
ciency of planetary systems formation to be in the range 




Cluster Size (N) 

Fig. 1. — Capture efBciency as a function of cluster size (Af*). The 
fraction of capture-formed systems is a decreasing function of the 
cluster size/mass, as expected from theoretical models; however, 
the dependence is weaker than the l/A^* dependence predic ted by 
models of soft binaries freeze-out ijMoeckel & Clarkell20lil ). The 
H and F models correspond to homogeneous and fractal models, 
respectively. 

of 0.036 < /pfP„^j X (/ffp/1) < 0.058, comparable to that 
of capture-formed stellar systems. 

Observations of the Taurus- Auriga star-forming region 
show a wide bina ry frequency of '-^ 5 — 10 % in the range 
1000 - 5000 AU ( Kraus et all 120 111 ), which is 2-3 times 
larger than found in the fractal models, and even higher 
compared with the homogeneous models. Note, how- 
ever, that for stars at these separations other processes 
beside dynamical capture are also likely to contribute to 
binary formation, and therefore direct comparison should 
be taken with caution. 

3.2. Properties of capture-formed systems 

In the following we discuss the distributions of various 
properties of the capture formed planetary and stellar 
systems. For simplicity we only present results from sim- 
ulations of Ni, = 100 clusters, which generally provide a 
good representation of the overall distributions we ob- 
serve. 

Mass function: Fig. [3] shows the fraction of 
capturcd-formed systems as function of the primary- 
component mass. As can be seen, the fraction is a mono- 
tonically increasing function of mass (which becomes sat- 
urated at high masses), for both stellar and planetary 
systems. 

SMA and closest approach distributions: The 

detailed distributions of SMAs and peri-center approach 
for various models are shown in Fig. (4] The orbits 
of captured stellar binaries and FFPs both follow quite 
similar distributions. The majority of the systems have 
SMAs in the typical range few x 10^ - 10'' AU, with thin 
tails extending towards lower and larger SMAs. Bina- 
ries with separations larger than ~ 10® AU are likely to 
be dissolved by the Galactic tide on timescales longer 
than our simulation, or through encounters with other 
field stars. Such systems are generally rare. The typical 
capture-formed planetary systems found in our homo- 
geneous cluster simulations follow an almost log-normal 
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TABLE 2 
Fractions of capture-formed systems 





Multiple 

Stars (%) 


Multiple 

Planets (%) 


Binary FFP 

(%) 


Binary star 
Planets (%) 


Triple star 
Planets (% ) 


^ Stellar/Planetary 


2/0 


3/0 


4/0 


1/1 


1/2 


1/3 


0/2 


2/1 


2/2 


3/1 


Model 






















Homogeneous 






















H_N10_P1 


10.3 


2.07 


0.12 


7.79 


0.71 


0.00 


0.01 


2.3 


0.18 


0.26 


H_N20_P1 


8.84 


1.49 


0.10 


6.51 


0.36 


0.01 


0.00 


1.6 


0.04 


0.12 


H_N50_P1 


6.59 


1.18 


0.10 


4.14 


0.18 


0.00 


0.00 


0.88 


0.03 


0.10 


H_N100_P1 


5.36 


0.86 


0.09 


3.01 


0.06 


0.01 


0.00 


0.56 


0.03 


0.08 


H_N200_P1 


4.24 


0.49 


0.03 


2.51 


0.07 


0.00 


0.00 


0.33 


0.01 


0.03 


H_N500_P1 


2.49 


0.24 


0.02 


1.61 


0.04 


0.00 


0.00 


0.14 


0.00 


0.02 


H_N1000_P1 


2.44 


0.16 


0.00 


1.41 


0.06 


0.00 


0.00 


0.14 


0.00 


0.01 
























H_N100_P1_50Mq 


5.43 


0.97 


0.07 


3.16 


0.14 


0.00 


0.00 


0.64 


0.02 


0.06 


H_N100_P0.5 


5.36 


1.15 


0.08 


1.53 


0.015 


0.00 


0.00 


0.35 


0.01 


0.04 


H_N100_P2 


4.88 


0.70 


0.03 


5.34 


0.31 


0.00 


0.01 


0.68 


0.05 


0.11 
























Fractal 






















F_N10_P1 


9.97 


2.65 


0.30 


8.71 


1.02 


0.04 


0.09 


3.69 


0.23 


0.55 


F_N20_P1 


8.69 


2.55 


0.35 


7.72 


0.68 


0.03 


0.02 


2.56 


0.16 


0.42 


F.N50.P1 


7.42 


1.90 


0.26 


6.06 


0.51 


0.01 


0.02 


1.57 


0.15 


0.35 


F.N100.P1 


6.79 


1.64 


0.17 


5.4 


0.37 


0.04 


0.01 


1.35 


0.09 


0.22 


F.N200.P1 


6.32 


0.91 


0.16 


5.0 


0.27 


0.00 


0.03 


0.93 


0.04 


0.12 


F.N500.P1 


5.24 


0.81 


0.12 


3.92 


0.22 


0.01 


0.02 


0.59 


0.04 


0.11 


F.N1000.P1 


5.4 


0.75 


0.12 


3.94 


0.18 


0.00 


0.02 


0.55 


0.03 


0.08 
























F.N100.P1.50Me 


6.52 


1.52 


0.18 


5.29 


0.32 


0.02 


0.03 


1.07 


0.07 


0.27 


F.N100.P0.5 


7.24 


1.71 


0.33 


2.78 


0.07 


0.01 


0.00 


0.72 


0.04 


0.1 


F.N100.P2 


5.87 


1.26 


0.23 


8.51 


0.79 


0.03 


0.10 


1.45 


0.17 


0.23 



distribution with a tail at low SMAs, and an upper cut- 
off mostly due to the Galactic tide, as mentioned above. 
The fractal cluster show a log-constant distribution at 
the mid range (lO^'^ — 10^'^ AU), i.e. a uniform distri- 
bution per logarithmic bin (Opic like distribution). Fig. 
|4|d shows the peri-center distribution, which extends as 
close as a few AU from the host stars (though these ex- 
treme cases are rare). Planets and stars captured into 
such close approaches can have a significant dynamical 
interactions with pre-existing planetary /protoplanetary 
systems even in dynamical timescales following their cap- 
ture. 

We note that there is a weak trend in the depen- 
dence of the SMA and peri-center distribution on cluster 
sizes. We find that closer binaries/planetary systems are 
formed in larger clusters (not shown); the overall distri- 
bution is generally well represented by the shown distri- 
bution in Fig. m 

Eccentricity distribution: We find the eccentricities 
of both planetary systems and stellar-binaries to follow 
a thermal distribution (see Fig. [2]), dN/de = 2e. 

Inclination distribution: In higher multiplicity sys- 
tems, we can explore the relative inclinations between 
inner and outer binaries of hierarchical systems. We find 
the inclinations distribution appears to be isotropic (uni- 
form in cos i). 

4. DISCUSSION 

The cluster dispersal scenario for the formation of wide 
orbit planetary systems is a natural extension for a sim- 
ilar model used to explain the origin of wide binaries. 
Our simulation results verify that FFFs can be captured 
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Fig. 2. — The cumulative eccentricity distribution of capture- 
formed planetary systems. The eccentricity distribution closely 
follows a thermal distribution (shown in solid line). 



through this process, with comparable efficiency. The 
capture process depends on both the cluster properties, 
as well as the masses of the specific stars. We find that 
the general properties of the wide planetary systems fol- 
low the same behavior as wide binar ies, which was ex- 
plored in detail in Papers I and in iMoeckel fc Clarki 
(|2011D . In the foUowing we focus on the specific im- 
plications for planetary systems. 
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Fig. 3. — The mass function of capture-formed system. Left: 
The mass function of the primary star /host in capture- formed bi- 
naries/planetary systems. Right: The fraction of capture-formed 
systems as a function of the primary mass. As can be seen, the cap- 
ture efficiency of both stellar and planetary companions increases 
with the host mass. The H and F models correspond to homoge- 
neous and fractal models, respectively. 



4.1. Planets and brown dwarfs at wide orbits 

The most basic result of our suggested FFP-capture 
scenario, verified in the detailed simulations, is the ro- 
bust production of wide-orbit planetary systems follow- 
ing cluster dispersal which include FFPs. The fraction 
of capture-formed wide planetary systems in our simu- 
lations is dependent on the their yet unknown original 
number fraction in the birth clusters. The recent mi- 
crolensing finding of FFPs, as well as theoretical predic- 
tions, suggest th at these fractions are high^ 0.5 — 2 FFPs 
per cluster star (jSumi et al.l 120111 ) . The FFP recapture 
mechanism could therefore play a major role in the pro- 
duction of wide-orbit planetary system. Such planets can 
regularly form around stars at small separation, then be 
ejected from the system, likely through scattering with 
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Fig. 4. — The semi-major axis and peri-center distribution func- 
tion of capture-formed stellar and planetary systems, in A'^^. = 100 
clusters. Top: The main figure shows the cumulative SMA dis- 
tribution, and the inset shows the differential SMA distribution. 
Capture formed systems from the the fractal clusters (F) are char- 
acterized by a log-uniform distribution (equal fraction in logarith- 
mic bins), where as the homogeneous clusters (H) give rise to a 
more log-normal distribution. Bottom: The same for closest ap- 
proach (peri-center) distributions. Note different axis scales. The 
distribution of different cluster sizes behave similarly with a gen- 
eral trend for the distribution in more massive clusters to peak at 
smaller separations. 

other planets, and then be recaptured by another clus- 
ter star following the birth-cluster dispersal. The rela- 
tion between such recaptured-planets and their new hosts 
could therefore be quite arbitrary. In particular the pair- 
ing of stellar and planetary properties such as metallicity, 
masses/sizes and spin-orbit relations should generally be 
random, or at most reflect correlations that may exist 
between the properties of any two objects in the birth 
cluster (e.g. related metallicity). Even if the latter case 
is true, the random pairing should be smeared out com- 
pared with the likely stronger correlation between pairs 
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of stars/planets originally formed in the same system. 

Planets could be captured into orbits of any eccentric- 
ity (through high eccentricities are preferred, through the 
thermal eccentricity distribution observed in our simula- 
tions), and a wide range of orbits, with typical SMAs 
extending over the range few x 10^ — lO'^ AU. At apo- 
center the most eccentric orbits could have twice as large 
separation as their SMA, and could therefore be more 
susceptible to disruption by the Galactic tide over time, 
if their SMA is comparabale to their Jacoby radius (Eq. 
[J). Systems with the smallest separation and/or the 
highest eccentricities, the peri-center approach, a(l — e), 
can become comparable to the separations of regularly 
formed planets at up to a few tens of AU. Objects on 
such orbits may strongly interact with any pre-existing 
planetary/protoplanetary system (see also next section). 

In-situ formation of planets (and particularly massive 
gas-giants) at wide orbits as found in our simulations, 
where the mass density in the protoplanetary disk be- 
come very small, is unlikely in the core- accretion model 
(jPodson- Robinson et al.ir2009f ). Ejection of planets from 
inner orbits to wide separations is possible. Such or- 
bits would initially be very eccentric and flyby encoun- 
ters with other stars would be required in order to put 
them into more circular orbits ( similar to the forma tion 
and evolution of the Oort cloud; iDuncan et al.1119871 ). In 
comparison, the capture scenario can naturally produce 
both low and high eccentricity orbits. Planets at wide 
orbits, however, are not likely to survive in a cluster en- 
vironment where close encounters with other stars can 
destroy the system (e.g. iParker fc Quanzir2012l and refer- 
ences therein). If such systems survive, or scattering oc- 
curs after the cluster dispersal it might be more difficult 
to currently distinguish between such a scenario and the 
suggested capture scenario introduced here. Neverthe- 
less, our models provide specific predictions for the dy- 
namical properties of such systems, and moreover, con- 
nect them to wide binary systems and population of free 
floating planets. The models therefore provide a wealth 
of potential signature s, which can distin guish them from 
different models (e.g. iVeras et al.ll2009l ). A smoking gun 
signature for the capture scenario may exist in the form 
of binary FFPs (see section 14. 4p . and the existence of 
planets in wide orbits around low mass stars, which are 
unlikely to form through in-situ planet formation sce- 
narios. In the gravita tional instability model for planet 
formation (|Bosj|2006[) massive gas planets may possibly 
form at wider orbits than in the core-accretion scenario, 
but it is not clear if the process could be efficient up 
to the very large separations as observed in our simula- 
tions. In any case, such processes are not likely to form 
low mass rocky planets at these separations. To sum- 
marize, in the FFP recapture model we expect planet of 
any mass to be captured into very wide orbits; systems 
which would be difficult to form in the currently sug- 
gested planet formation models. We also note that the 
fraction of capture-formed planetary systems, as well as 
the planetary masses distribution reficct the FFP frac- 
tion and mass distribution in the birth-clusters, and are 
likely to show similar orbital properties as wide stellar 
binaries. 

4.2. Multiple planets at wide orbits 



In principle more than a single planet can be cap- 
tured into the same planetary system, through indepen- 
dent captures. The fraction of such multiple captured- 
planet systems is therefore expected to be proportional 
to the probability for multiple independent captures, i.e. 
/2p = /p for capture of two planets. Table 2 show the 
fraction of planetary systems with n = 1, 2 and 3, which 
is generally consistent with the expectation. Capture of 
high (> 2) multiplicity stellar systems is expected to fol- 
low the same behavior, however, capture of an additional 
star also contributes to the total mass of the system, and 
therefore the likelihood of a consecutive capture is some- 
what increased. This is verified in the higher fraction 
of high multiplicity stellar systems in comparison to the 
planetary systems (when normalized to the fraction of 
binary/single planet systems). 

Though this scenario can explain even multiple planet 
systems at wide orbits, it is unlikely to form a large num- 
ber of such systems, given the rapidly decreasing prob- 
ability for such occurrence. It is also unlikely to form 
planetary systems in mean motion resonance, which re- 
quire a continuous migration of planetary orbits through 
phase space; these are likely to form through migration 
in disks. 

Naturally, stars hosting primordial in-situ formed plan- 
etary systems or protoplanetary disk can capture addi- 
tional planets following the cluster dispersal. Given the 
random orientation of the capture processes, these ad- 
ditional 'guest' captured planets are likely to be mis- 
aligned with any existing protoplanetary disk or pre- 
existing planetary system. In some cases, the capture 
of an additional object into a wide orbit may lead to the 
perturbation or even the destabilization of an existing 
planetary system. In particular, an object captured at 
high inclination may induce long term secular evolution 
of the inner pre- existing system (e.g. th rough Kozai- 
Lidov evolution; lKozailll963 lLidovlTl962l ) . and poten- 
tially even destabilize them (though the latter possibility 
is far more likely in the case of a stellar capture). 

4.3. Planets in binary systems: circumstellar and 
circumbinary orbits 

Primordial binaries existing before the cluster dispersal 
could themselves capture additional stellar or planetary 
objects, forming high multiplicity stellar systems of bi- 
nary planets. This mechanism is not much different than 
the capture of planet around single stars (aside from the 
potentially much higher mass of the system, in the case 
of binary stars, which induces a higher capture probabil- 
ity, see Fig. Nevertheless, three body systems could 
become unstable, and therefore the capture of additional 
objects by a binary may limit the range of orbital phase- 
space available, at least in those cases where the initial 
binary is wide. The current simulations presented here 
did not include primordial binaries, which would be a 
natural extension of our models. We expect such sys- 
tems to generally have a higher capture efficiency due to 
their higher mass. 

Binary planets could also form through a double cap- 
ture; either a planet capture followed by a stellar cap- 
ture or a stellar capture, forming a binary, followed by 
a planetary capture. The planetary capture in the lat- 
ter case, follows a similar evolution as planetary capture 
around primordial binaries; however, the inner stellar bi- 
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nary would be typically much wider than primordial bi- 
naries in the birth cluster. The former process where a 
planet is captured first, may lead to a circumstellar bi- 
nary planet, as the stellar binary could be captured at 
a much wider orbit, not perturbing the inner planetary 
orbit. In this case, circumbinary planetary systems is 
less likely to form, as a close single-binary encounter of 
a star within the the planetary orbit (i.e. where the bi- 
nary is the star-planet system), would more likely lead 
to its ejection process (typically the lowest mass objects 
are ejected in such encounters). Table 2 shows the frac- 
tion of binary planets (combining both circumstellar and 
circumbinary orbits) formed in our simulations. 

4.4. Free floating planet-binaries 

Larger stellar masses allow for binding of companions 
with higher relative velocities, as reflected in the higher 
capture rate by more massive stars. Nevertheless, planets 
could capture companion planets without any host stars. 
These free floating double planets or FFP-binaries, are 
generally rare (see Table 2). The formation of such sys- 
tems (or similarly planet capture around brown dwarfs) 
through typical planet formation scenarios rather than 
through dynamical capture processes is unlikely (cur- 
rent models require the existence of a host stars around 
which planet forms in a massive enough protoplanetary 
disk; see also di scussion on planet and star formation 
at wi de orbits in IBeiar et"al]|2008t IKraus fc Hillenbrandl 
|2009() . Therefore, the finding of even a single such sys- 
tem could serve as a potential smoking gun signature of 
the planetary capture scenario discussed here. Given the 
increased capture probability with system mass, and the 
better detectability of more massive planets, the most 
likely FFP-binaries to be observed would lie in the higher 
end of the planetary mass reg ime. Most interestingly, 
IJavawardhana fc Ivanovl ()2006f ) discovered such a sys- 
tems of two planetary mass companions of a fewxlOMj 
at wide, 240 AU separation, followed by the discovery of 
similar even wider syst ems by other groups (|Beiar et al.l 
[200llBiller eraI|[20lTI ) . Our suggested scenario can pro- 
vide a natural explanation for these otherwise puzzling 
systems. 

4.5. Wide orbit planets around compact objects 

The timescales for cluster dispersal could be relatively 
short, and much shorter than the main-sequence lifetimes 
of most stars. Nevertheless, the lifetime of stars more 
massive that 5 Mq could be comparable to the dis- 
persal time scales of the host cluster. These stars could 
therefore become compact objects before the cluster dis- 
perses, and can later capture planets like any other clus- 
ter star. Note, however that this is likely to be true only 
for black holes and white dwarfs; the natal kicks of neu- 
tron stars are likely to ejects them from their host cluster 
upon their formation. The mass loss of stars prior to the 
formation of the compact object can eject pre-existing 
planets at wide orbits (e.g. Veras et al. 2012, submitted); 
especially in the case of the prompt mass loss in a super- 
nova. The capture channel suggest that even such objects 
could host wide-orbit planetary systems. Interestingly, 
the dispersal timescale of a cluster defines a lower limit 
on the mass of WD progenitors that could acquire plan- 
etary/binary companions through capture. Lower mass 



WDs with such wide orbits companion would have had 
to capture them while still on the MS, and would have a 
higher chance of losing these companions as they evolve 
and lose mass. The age cut-off may therefore be reflected 
by a lower fraction of low mass WDs with wide orbit 
stellar or planetary companions. Masses of stellar black 
holes could typically be in the range of 5 — 15 Mq; given 
the higher probability for planetary and stellar capture 
for massive stars (see Fig. [3]), we expect more than half 
of all stellar black holes to have a wide orbit captured 
stellar companion, and more than 5 — 10 x (/ffp/1) % 
of them to host a planetary companion (our additional 
simulations of clusters containing stars as massive as 20 
Mq , not shown here, suggest that the capture probability 
for stellar companions could approach unity for rrii, > 10 
Mq). Massive white dwarfs (mwD ^ 1 Mq) are expected 
to capture ~ 0.2 less companions compared with stellar 
black holes. 

5. SUMMARY 

In this paper we presented a novel channel for the for- 
mation of planetary systems at very wide orbits (typ- 
ically ~ few X 10^ — 10^ AU), through the dynamical 
recapture of free floating planets during a cluster disper- 
sal. This mechanism, originally suggested as a route for 
the formation of wide orbit binaries, was extended to the 
planetary regime. We simulated dispersing clusters with 
population of FFPs, and followed the formation of plan- 
etary systems in the simulations. We find that the effi- 
ciency of capture-formation of planetary systems ranges 
between 1 — 9 % depending on cluster size and struc- 
ture. We predict and overall fraction of 3 — 6 x (/ffp/1) 
% of all stars (accounting for the birth cluster size dis- 
tribution) to host a wide orbit planetary companion at 
separations > 100 AU, comparable to the expected frac- 
tion of wide stellar binaries. The orbital properties of 
capture-formed planetary systems resemble those of wide 
binaries. The fraction of capture-formed planetary sys- 
tems (as well as capture-formed binaries) increases with 
the host mass. Future searches for wide-orbit planetary 
companions should therefore favor more massive stars. 
Although captures are biased towards higher mass hosts, 
planets can be captured around hosts of any mass (with 
stars like our Sun having a 3 — 5 % probability to have 
captured a wide orbit planet). In fact, wc find that, al- 
though rare, planets could even be captured around other 
planets (or brown dwarfs) to form free fioating planet- 
binaries. In addition, planets can be captured around 
compact objects, if these formed prior to the cluster dis- 
solution; we find the stellar black holes may have a high 
probability of capturing stellar and planetary compan- 
ions in wide orbits. Finally, we note that the properties 
of captured planets should not correlate with their host 
star; single and multiple planets could be captured into 
arbitrary orbits, generally non-coplanar with each other 
(or with the stellar host spin, or pre-existing planetary 
system). Nevertheless, some secondary planet- host star 
metallicity correlation may still exist due to the common 
origin of the stars from same birth cluster. 
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